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Abstract 

We present a supersymmetric model of quark and lepton based on S'4 x Z3 x 
flavor symmetry. The 5*4 symmetry is broken down to Klein four and Z3 subgroups 
in the neutrino and the charged lepton sectors respectively. Tri-Bimaximal mixing and 
the charged lepton mass hierarchies are reproduced simultaneously at leading order. 
Moreover, a realistic pattern of quark masses and mixing angles is generated with the 
exception of the mixing angle between the first two generations, which requires a small 
accidental enhancement. It is remarkable that the mass hierarchies are controlled by 
the spontaneous breaking of flavor symmetry in our model. The next to leading order 
contributions are studied, all the fermion masses and mixing angles receive corrections 
of relative order with respect to the leading order results. The phenomenological 
consequences of the model are analyzed, the neutrino mass spectrum can be normal 
hierarchy or inverted hierarchy, and the combined measurement of the Oz^2/3 decay ef- 
fective mass mpp and the lightest neutrino mass can distinguish the normal hierarchy 
from the inverted hierarchy. 



^e-mail address: dinggj@ustc.edu.cn 



1 Introduction 



Neutrino has provided us a good window to the new physics beyond the Standard ModeL 
Neutrino oscillation experiments have provided solid evidence that neutrinos have small but 
non-zero masses. Global data fit to the current neutrino oscillation data demonstrates that 
the mixing pattern in the leptonic sector is so different from the one in the quark sector. Two 
independent fits for the mixing angles and the mass squared differences are listed in Table [H 



parameter 


Ref. [1] 
best fitilcr 3cr interval 


Ref. [2] 

best fitilcr 3a interval 


Amii [lO-^eV^] 
Ami^ [lO-^eV^] 

sin^ 9i2 
sin^ 623 
sin^ 013 


7.65l[J:| 7.05-8.34 

±2AOto\l ±(2.07-2.75) 

0.304l°:°^2 0.25-0.37 
O.dOtofe 0.36-0.67 
0.01^;^:j^i? < 0.056 


7.67^[J|? 7.07-8.34 

-2.39 ±0.12 -(2.02-2.79) 
+2.49 ±0.12 +(2.13-2.88) 
0.32ll°:°23 0.26-0.40 

0.47t[5;J]^ 0.33-0.64 

0.003 + 0.015 < 0.049 



Table 1: Three flavour neutrino oscillation parameters from two global data fits [1,2]. 



As is obvious, the current neutrino oscillation data is remarkably compatible with the so 
called Tri-Bimaximal (TB) mixing pattern [3], which suggests the following mixing pattern 
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sin"" 6i2,TB = 7;, sin^ 023,tb = -, sin^ 6*13,^5 = 



These values lie in the la range of global data analysis shown in Table W- Correspondingly, 
the leptonic Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is given by 
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where 021 and 031 are the Majorana CP violating phases, and Utb is given by 
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The mixing in the quark sector is described by the famous CKM matrix [4], and there is 
large mass hierarchies within the quarks and charged leptons sectors respectively [5]. The 
origin of the observed fermion mass hierarchies and flavor mixings is a great puzzle in particle 
physics. Nowadays promising candidates for understanding such issue are the models based on 
spontaneously breaking flavor symmetry, various models based on discrete or continuous flavor 
symmetry have been proposed so far [6,7]. Recently it was found that flavor symmetry based 
on discrete group is particularly suitable to reproduce specific mixing pattern at leading order 



^svY? ByiTB is exactly within the \a range of the second global data fit, whereas it slightly above the \a 
up limit of the first fit. 
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[8] . The A4 models are especially attractive, it has received considerable interest in the recent 
past [9-14]. So far various flavor models have been proposed, and their phenomenological 
consequences were analyzed [15-18]. These models assumed that A4 symmetry is realized at 
a high energy scale, the lepton fields transform nontrivially under the symmetry group, and 
the flavor symmetry is spontaneously broken by a set of flavons with the vacuum expectation 
value (VEV) along a specific direction. The misalignment in the flavor space between the 
charged lepton and the neutrino sectors results in the TB lepton mixing. 

If extend the A4 symmetry to the quark sector, the quark mixing matrix Vckm turns 
out to be unity matrix at leading order [11], However, the subleading contributions of the 
higher dimensional operators are too small to provide large enough deviations of Vckm from 
the identity matrix. The possible ways of resolving this issue are to consider new sources 
of symmetry breaking or enlarge the symmetry group. Two discrete groups T' [19-26] and 
5*4 [27-32] are found to be promising, both groups have two dimensional irreducible represen- 
tation, which is very useful to describing the quark sector. The S4 symmetry is particularly 
interesting, 5*4 as a horizontal symmetry group has been proposed long ago [33], and some 
models with different purposes have been built [34]. Recently it was claimed to be minimal 
flavor group capable of yielding the TB mixing without flne tuning [35-37]. However, Grimus 
et al. were against this point [38]. 

In this work, we build a SUSY model based on S4 x x Z4 flavor group, the neutrino 
mass is generated via the conventional type I See-Saw mechanism [39]. Our model naturally 
produces the TB mixing and the charged lepton mass hierarchy at leading order. Furthermore, 
we extend the model to the quark sector, the realistic patterns of quark masses and mixing 
angles are generated. In our model the mass hierarchies are controlled by the spontaneous 
breaking of the flavor symmetry instead of the Froggatt-Nielsen (FN) mechanism [40]. 

This article is organized as follows. Section 2 is the group theory of 6*4 group, where 
the subgroup, the equivalent class, and the representation of S4 are presented. In section 3 
we justify the vacuum alignment of our model in the supersymmetric limit. In section 4 we 
present our model in both the lepton and quark sectors, its basic features and theoretical 
predictions are discussed. In section 5 we analyze the phenomenological implications of the 
model in detail, which include the mass spectrum, neutrinoless double beta decay and the 
Majorana CP violating phases etc. The corrections induced by the next to leading order 
terms are studied in section 6. Finally we summarize our results in the conclusion section. 



2 The discrete group S4 

Si is the permutation group of 4 objects. The group has 24 distinct elements, and it can 
be generated by two elements S and T obeying the relations 

S^ = T^ = 1, ST^S = T (4) 

Without loss of generality, we could choose 

S = (1234), T = (123) (5) 

where the cycle (1234) denotes the permutation (1,2,3,4) — (2,3,4,1), and (123) means 
(1,2,3,4) — > (2,3,1,4). The 24 elements belong to 5 conjugate classes and are generated 
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from S and T as follows 



Ci:l 

C2 : STS^ = (12), TSTS'^ = (13), ST^ = (14), ^^^^ = (23), TST = (24), T^^ = (34) 
C3 : TS^T^ = (12)(34), = (13)(24), T^S^T = (14)(23) 

C4 : T = (123), = (132), T^S^ = (124), S^T = (142), S^TS^ = (134), STS = (143), 

^2^2 ^ (^234), = (243) 
€5:8 = (1234), T^^T = (1243), ST = (1324), T5 = (1342), TST^ = (1423), = (1432) 

The structure of the group 5*4 is rather rich, it has thirty proper subgroups of orders 1, 2, 3, 
4, 6, 8, 12 or 24. Concretely, the subgroups of S4 cirG clS follows 

1. The trivial group only consisting of the unit element. 

2. Six two-element subgroups generated by a transposition of the form {1, (ij)} with i ^ j 

= {1,5T52}, = {1,TSTS^}, = {l,ST^}, = {l,S^TS}, Hf^ = 
{1,T5T} and Hf^ = {l,T^S} 

3. Three two-element subgroups generated by a double transition of the form {1, {ij){kl)} 
with i j ^ k ^ I 

hP = {1,TS^T^}, hP = {1,S^}, hP = {1,T^S^T} 

4. Four subgroups of order three, which is spanned by a three-cycle 

hP = {1,T,T^},hP = {1,T^S^S^T},hP = {1,S^TS^STS},hP = {1,S^T^,TS^} 

5. The four-element subgroups generated by a four-cycle, they are of the form {1, g, g"^, g^} 
with g any four-cycle 

hP = {1,S,S^,S^}, hP = {1,TS,T^ST,T^S'^T},hP = {1,TST'^,ST,TS'^T^} 

6. The four-element subgroups generated by two disjoint transpositions, which is isomor- 
phic to Klein four group 

hP = {1, STS"^, T^S, TS^T^}, hP = {1, TSTS"^, TST, S^}, Hf = {1, ST^, S^TS, T^S^T} 

7. The order four subgroup comprising of the identity and three double transitions, which 
is isomorphic to Klein four group 

hP = {l,TS^T\S\T^S^T} 

8. Four subgroups of order six, which is isomorphic to 53. They are the permutation 
groups of any three of the four objects, leaving the fourth invariant 

hP = {1,STS^,TSTS^,S^TS,T,T^}, = {1, STS^, ST^,TST,T^S\ S^T}, 
hP = {1,TSTS^ ST^, T^S, S^TS\ STS}, = {1, S^TS, TST, T^S, S^T^, TS^} 
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C2 
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C4 
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-1 
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-1 
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-1 





3i 


3 
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-1 





-1 


32 


3 


-1 


-1 
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Table 2: Character table of the ^4 group. 12^ denotes the number of the elements contained in the class Ci, 
and is the order of the elements of Ci . 



9. Three eight-element subgroups, which is isomorphic to D4 

H^^'' = {1,TSTS^ TST, S, S\ TS^T^, S^, T^S^T], 
i/f ^ = {1, STS^, T^S, ST, TST^, TS^T^, S^, T^S^T} 
H^^^ = {1, ST^, S^TS, T^ST, TS, TS^T^, S\ T^S^T} 

10. The alternating group A4 

= {1, TS^T^, S"^, T^S^T, T, T2, T^S"^, S^T, S^TS^, STS, S^T^, TS^} 

11. The whole group 

In particular, Hi and A4 are the invariant subgroups of 5*4. Since the number of the un- 
equivalent irreducible representation is equal to the number of class, the 6*4 group has five 
irreducible representations: li, I2, 2, 3i and 32- li is the identity representation and I2 is the 
antisymmetric one. The Young diagram for the two dimensional representations is self associ- 
ated, and the Young diagrams corresponding to the three dimensional representations 3i and 
32 are associated Young diagrams. For the same group element, the representation matrices 
of 3i and 32 are exactly the same if the element is an even permutation. Whereas the overall 
signs are opposite if the group element is an odd permutation. It is notable that S4 together 
with T' is the smallest group containing one, two and three dimensional representations. The 
character table of S'4 group is shown in Table [21 

From the character table of the S'4 group, we can straightforwardly obtain the multiplica- 
tion rules between the various representations 

Ij Ij = l((i+j) mod 2)+l, Ij ® 2 = 2, Ij ® 3j = 3((j+j) mod 2)+l 

2®2 = li©l2©2, 2®3i = 3i©32, 3^ ® 3^ = li © 2 © 3i © 32, 

3i©32 = l2©2©32©32, with «,j = 1,2 (6) 

The explicit representation matrices of the generators S, T and other group elements for the 
five irreducible representations are listed in Appendix A. From these representation matrices, 
one can explicitly calculate the Clebsch-Gordan coefficients for the decomposition of the 
product representations, and the same results as those in Ref. [29] are obtained. 
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Table 3: The transformation rules of the matter fields and the flavons under the symmetry groups 6*4, Z3 and 
Z4. Lo is the third root of unity, i.e. w = e*T^ = (—1 + iV3)/2. We denote Ql = {Qi, Q2Y which are doublets 
of 54, where Qi = {u,dY and Q2 = (c, s)* are the electroweak SU(2) doublets of the first two generations. 
(5.3 = (t, hy is the electroweak SU(2) doublet of the third generation. 

3 Field content and the vacuum alignment 

The model is supersymmetric and based on the discrete symmetry 5*4 x Z3 x Z^. Super- 
symmetry is introduced in order to simphfy the discussion of the vacuum ahgnment. The 
component controls the mixing angles, the auxiliary Z3 symmetry guarantees the misalign- 
ment in flavor space between the neutrino and the charged lepton mass eigenstates, and the 

component is crucial to eliminating the unwanted couplings and reproducing the observed 
mass hierarchy. The fields of the model and their classification under the flavor symmetry are 
shown in Table [3l where two Higgses doublets hu,d of the minimal supersymmetric standard 
model are present. If the 6*4 flavor symmetry is preserved until the electroweak scale, then 
all the fermions would be massless. Therefore 5*4 symmetry should be broken by the suitable 
flavon fields, which are standard model singlets. Another critical issue of the flavor model 
building is the vacuum alignment, a global continuous symmetry is exploited to sim- 

plify the vacuum alignment problem. This symmetry is broken to the discrete R parity once 
we include the gaugino mass in the model. The matter fields carry -|-1 R-charge, the Higgses 
and the flavon supermultiplets have R-charge 0. The spontaneous breaking of 6*4 symmetry 
can be implemented by introducing a new set of multiplets, the driving fields carrying 2 unit 
R-charge. Consequently the driving fields enter linearly into the superpotential. The suitable 
driving fields and their transformation properties are shown in Table HI In the following, we 
will discuss the minimization of the scalar potential in the supersymmetric limit. At the lead- 
ing order, the most general superpotential dependent on the driving fields, which is invariant 
under the flavor symmetry group 6*4 x Z3 x Z4, is given by 



Fields 


^0 




9^ 




00 


AO 


S4 


3i 


I2 


ll 


2 


32 


ll 


Z3 


1 


1 


1 








Z4 


-1 


-1 


1 


1 


1 


1 



Table 4: The driving fields and their transformation properties under the flavor group ^4 x Z3 x Z4. 



h{ri\m)2)u + /2(r/°(00)2)i. + /3(0°(#)3ju + /^iA°A2 + h^^^^iii),, + h,^\<P<P),, (7) 
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where the subscript li denotes the contraction in li, similar rule applies to other subscripts 
I2, 2, 3i and 82. In the SUSY limit, the vacuum configuration is determined by the vanishing 
of the derivative of w„ with respect to each component of the driving fields 



'^giivi - + 2fi'2(x? - X2X3) + 93(^2X3 - V3X2) = 

25'i('/^2 - V^iV^s) + 25f2(X2 - X1X3) + 5'3('^3Xi - V'lXs) = 

"^gii^l - V1V2) + 2fi'2(x3 - X1X2) + g3{.viX2 - V2X1) = 

9a{^iXi + ^2X3 + ^3X2) = 
M| + Ke^ = 



89^ 

This set of equations admit the solution 

((^) = (0, tv, 0), (x) = (0, v^, 0), {9) = ve (9) 

with 

2 92 2 2 ^'^0 / \ 

v^p = v^, Vg = , undetermined (10) 

9i 1^ 

From the driving superpotential w^, we can also derive the equations from which to extract 
the vacuum expectation values of rj, and A 



drfl 

Wi 

W2 
dWy 

dWy 



/l^l' + /2(03 + 20102) =0 
/ir/2' + /2(0^ + 20103) =0 
/3(?7l02 - ?7203) = 
/3('7l</'l - V2(l)2) = 

f3ivih - V2<i)i) = 

/iiA^ + 2Mir/2 + /i3(0? + 20203) = (11) 



9A0 

The solution to the above six equations is 

iv) = {Vr„Vrj), (0) = {V^,V^,V^), (A) = Va (12) 

with the conditions 

2 2 2 fl^3 — 2 /'2^2 2 / \ 

= "TTT'^r?' "^A = TT "^jyi '^v uudetcrmmed (13) 



3/2 ^ hh 



The vacuum expectation values (VEVs) of the fiavons can be very large, much larger than the 
electroweak scale, and we expect that all the VEVs are of a common order of magnitude. This 
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is a very common assumption in the flavor model building, which guarantees the reasonability 
of the subsequent perturbative expansion in inverse power of the cutoff scale A. Acting on 
the vacuum configurations of Eq.Q and Eq. fll2l) with the elements of the flavor symmetry 
group 5*4, we can see that the VEVs of t] and (p are invariant under four elements 1, TST, 
TSTS"^ and 5*^, which exactly constitute the Klein four group On the contrary, the 

VEVs of ip and x break 5*4 completely. Under the action of T or T^, the directions of ((/?) 
and (x) are invariant except an overall phase. Considering the enlarged group 5*4 x Z3, the 
vacuum configuration Eq.(9) preserves the subgroup generated by uT, which is defined as 
the simultaneous transformation of T G 6*4 and a; G Z3. As we shall see later that the 5*4 
flavor symmetry is spontaneously broken down by the VEVs of rj and (p in the neutrino sector 
at the leading order (LO), and it is broken down by the VEVs of ip and x iii the charged lepton 

sector. Whereas both rj, and if, x are involved in generating the quark masses. The 5*4 flavor 

(5) 

symmetry is broken into the Klein four symmetry and the symmetry generated by 
T in the neutrino and the charged lepton sector respectively at LO. This symmetry breaking 
chain is crucial to generating the TB mixing. 



4 The model with S4 x x Z4 flavor symmetry 

In this section we shall propose a concise supersymmetric (SUSY) model based on 5*4 x 
Z^ X Z4 flavor symmetry with the vacuum alignment of Eq.(9) and Eq.(12). 



4.1 Charged leptons 

The charged lepton masses are described by the following superpotential 

^^ = ^ e^(^<^)ii(w)ii^d + ^ e'((£(^)2(<^<^)2)ii/id + ^ e^{{£<f)3,{<f<f)3,)uhd 
+ ^ e%{iipUxx)2)i^d + ^ e%iiphAxx)3juhd + ^ e^{{ix)2{^Muhd 

In the above superpotential W£, for each charged lepton, only the lowest order operators in 
the expansion in powers of 1/A are displayed explicitly. Dots stand for higher dimensional 
operators. Note that the auxiliary Z4 symmetry imposes different powers of (p and x for the 
electron, mu and tau terms. At LO only the tau mass is generated, the muon and the electron 
masses are generated by high order contributions. After the flavor symmetry breaking and 
the electroweak symmetry breaking, the charged leptons acquire masses, and wi becomes 



We 



{ye2 - 22/e3)T^ + (-2/e4 + 22/e5)T^ + iVer " '^Ves)^ + (-^/eQ + 2?/, 



^3 i V ye^ i -»eo/^3 > \ye, -yeoj ^3 > v yea i -yelOj ^3 



Vde e 



3 

Vej^ Vde e + y^,^;^ Vd^x /i + Vrj^ v^t t (15) 
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where Vd = {hd), Ve = ye2 - 2?/e3 + (-l/e4 + 2l/e5) J + {Vel " '^Ves)^ + (-|/e9 + 2|/el0 

= 2y'^. As a result, the charged lepton mass matrix is diagonal at LO 



rrii 



( 



\ 



Ve ^3 



















\ 



J 



Vd 



and 



(16) 



It is obvious that the hermitian matrix m\mi is invariant under both T and T displayed in 
the Appendix A, i.e., 

T'^m\m^T = mlmi (17) 

Conversely, the general matrix invariant under T and must be diagonal. Consequently 
the 5*4 symmetry is broken to the subgroup h[^^ = in the charged lepton sector. The 
charged lepton masses can be read out directly as 

,3 



VeJ^Vd 



m. 



A2 



Vd 



Vu, 
Vr-^Vd 



we notice that the charged lepton mass hierarchies are naturally generated by the spontaneous 
symmetry breaking of 5*4 symmetry without exploiting the FN mechanism [11]. Using the 
experimental data on the ratio of the lepton masses, one can estimate the order of magnitude 
of v^/A and v^/A. Assuming that the coefficients ye, and yr are of 0{1), we obtain 



me V 

— ~ — ^ ~ 3 X 10 

m-r A'^ 



— - ~ ~ 6 X 10 
rrir A 



-2 



(19) 



Obviously the solution to the above equations is 



^) ~ (±1.73 X 10-2, 6 X 10-2) 
' A A ^ ^ / 



(20) 



we see that the amplitudes of both f A and v^/ A are roughly of the same order about 0{\ 
where Ac is the Cabibbo angle. 



2) 



4.2 Neutrinos 

The superpotential contributing to the neutrino mass is as follows 

= ^-f{{^'^)2V)uhu + ^-f{{v^^hA)uhu + \m{v'v%, + ... (21) 

where dots denote the higher order contributions, M is a constant with dimension of mass, 
and the factor | is a normalization factor for convenience. The first two terms in Eq.( |2T]) 
determine the neutrino Dirac mass matrix, and the third term is Majorana mass term. After 
electroweak and 5*4 symmetry breaking, we obtain the following LO contributions to the 
neutrino Dirac and Majorana mass matrices 

I 2h a-b a-b \ / M \ 

m^=la-ba + 2b -b \ Vu, Mn = I M (22) 
\ a-b -b a + 2b I \0M0/ 
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where Vu = {hu), a = Uvi^ and b = yu2-^- We notice that the Dirac mass matrix is symmetric 
and it is controlled by two parameters a and b. The eigenvalues of the Majorana matrix 
are given by 

Ml = M, Ms = M, Ms = -M (23) 

The right handed neutrino masses are exactly degenerate, this is a remarkable feature of our 
model. Integrating out the heavy degrees of freedom, we get the light neutrino mass matrix, 
which is given by the famous See-Saw relation 



2a2 + - 4a6 - Sfe^ + 2ab - 36^ + lab 



m„ = -(m^fM^^m^ = --^\ - 3b^ + 2ab - 3b^ - 4ab 2a^ + 6b^ + 2ab 

a^-3b^ + 2ab 2a^ + Qb^ + 2ab - 3b^ - Aab 

(24) 

The above light neutrino mass matrix m,^ is 2 3 invariant and it satisfies the magic 
symmetry (m;^)ii + (m;^)i3 = {171^)22 + {^11^)23. Therefore it is exactly diagonalized by the TB 
mixing 

U'^m^Uy = diag(mi, m2, mg) (25) 
The unitary matrix U^, is written as 

= Utb diag(e-*"^/2^ e'^^^/^^ g"'"^/') (26) 

The phases ai, 02 and are given by 

ai = arg(-(a-36)VM) 
a2 = arg{-Aa^/M) 

as = arg((a + 36)VM) (27) 
nil, ^2 and 1713 in Eq.(2[S]) are the light neutrino masses. 



nil = \ {a — 3by 
m2=4|a2' " 



|M| 



|M| 

mg = |(a + 36)^1^ (2J 



\M\ 

Concerning the neutrinos, the 6*4 symmetry is spontaneously broken by the VEVs of 77 
and (p at the LO. since both (77) and (0) are invariant under the actions of TSTS"^, TST 
and 5*2, the flavor symmetry 5*4 is broken down to the Klein four subgroup = H^'^ = 
{1,TSTS^,TST, S"^} in the neutrino sector. We can straightforwardly check that the light 
neutrino mass matrix mj, is really invariant under TSTS"^, TST and 5*^. On the contrary, 
the most general neutrino mass matrix invariant under the Klein four group is given by 

mil rni2 nii2 

my = I mi2 m22 mn + m^ - m22 \ (29) 

mi2 mil + mi2 - m22 ^22 
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where mn, mu and 77122 are arbitrary parameters. In the present model, the hght neutrino 
mass matrix is given by Eg. fl24l) . which is a particular version of the neutrino mass matrix in 
Eq. fl29p . Since only two parameters a and b are involved in our model, additional constraint 
has to be satisfied, i.e. 3m^^ + Amumn — 4m22?7iii — Smfg — 4m22 = 0, which is generally not 
implied by the invariance under G^. This is because that in our model the fields which break 
5*4 are a doublet rj and a triplet 0, there are no further fiavons transforming as li or 82 which 
couple to the neutrino sector. 

In short summary, at the LO the S4 fiavor symmetry is broken down to and Klein four 
subgroup in the charged lepton and neutrino sector respectively. We have obtained a diagonal 
and hierarchical charged lepton mass matrix, the heavy neutrino masses are degenerate, and 
the neutrino mixing matrix is exactly the TB matrix. 



4.3 Effective operators 

In the previous section, the neutrinos acquire masses via the See-Saw mechanism. It is 
interesting to note that higher dimension Weinberg operator cloud also contribute to the 
neutrino mass directly, which may correspond to exchanging some heavy particles rather 
than the right handed neutrinos u'^. In the present model, these effective light neutrino mass 
operators are 

+^mjhu)2{<l>')2)u + ^mjK)M%,)^, + ^{{ihJKh,{v<l>K)u (30) 

With the vacuum configurations displayed in Eq. fll2l) . the high dimension operators wf,^^ leads 
to the following effective light neutrino mass matrix 



a + 27 P — ■y P — 1 

<"= I P-1 /5 + 27 «-7 ) ^ (31) 
/3 — 7 a — 7 f3 + 2y 



where 



" = 2a;^ + 42/1^ + 62:1-^ 
/3 = 2|/2^ + 6.2| 

7 = 4^ (32) 

Obviously rrfj^ have the same texture as that in Eq.(29), and it is remarkable that this mass 
matrix is diagonalized by TB matrix, 

Utb'^I^^^tb = diag{m'l;^-^ ,'m'2-^^ ,'m':/^) (33) 

where rrfi^ , ni^^ and ni^^ are the effective light neutrino masses coming from the above 
high dimension Weinberg operators, they are given by 

mf^ = (a-/5 + 37)^ 
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mf^ = (a + 2/5)^ 

m^^^ = + + 37) I (34) 

If we consider the parameters yui^v2 ~ C^(l)) yi,2 ~ C^(l)) ^1,2,3 ~ C'(l) and x ~ w; ~ C'(l), 
we get the ratio 

rn/^ M , , 

— r 35 

We see that the importance of Weinberg operators depends on the relative size of M and 
A. Since we have assumed that the hght neutrino masses mainly come from the See-Saw 
mechanism, the right handed neutrino mass M should be much smaller than the cutoff scale 
A. In the context of a grand unified theory, this corresponds to the requirement that M is 
of order 0{Mgut) rather than of 0{Mpianck)- In some flavor models, right handed neutrino 
masses are required to be below the cutoff A as well, in order to reproduce the experimental 
value of the small parameter Am^^j/Am^^^ [14,32]. Another convenient way of suppressing 
the contributions of the effective operators is to introduce auxiliary symmetry further, so that 
the Weinberg operators arise at much higher order and its contributions can be neglected. 

4.4 Extension to the quark sector 

The Yukawa superpotentials in the quark sector are 

Wq = Wu + Wd (36) 

In the up quark sector, we have 

3 2 , 3 

j=l i=l i=l 

2 , 8 

i=l i=l 



2 



+ Y.-§u^Qz{0^'\.K + - (37) 

i=l 



where 



(9(1) = {^^^^x.xx] 

0(3) = x4>\ <^#^ X^V, V'^A^, X0A2} 

(9(^) = {^^^xn (38) 

The superpotentials contributing to the down quark masses are as follows 

i=l i=l 1=1 
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i=l 



1=1 



A3 



i=l 



A3 



A3 



i=l 



A4 



(39) 



where 



(9(5) 
0(6) 
0(7) 



W^V, V^V, X^^, xV, V^X^, 



(40) 



Since the quantum numbers of 6^ and are exactly the same, as is obvious from Table El 
there are no fundamental distinctions between 6^ and s'^, we have defined b'^ as the one which 
couples to QsOhd in the superpotential Wd- We notice that both the supermultiplets (p, x and 
r), (f), which control the flavor symmetry breaking in the charged lepton and neutrino sectors 
respectively, couple to the quarks. Consequently the 6*4 flavor symmetry is completely broken 
in the quark sector. By recalling the vacuum configuration in Eq.(l9]) and Eq.( |T2i) . we can 
write down the mass matrices for the up and down quarks 





yii A4 


(n) -o^^l 

yi2 A" 


,,(") < 
yi3 A3 


\ 


rriu = 


(n) ^-^A 

y2i A3 


y22 A2 


y23 A 






\ 2/31 A^ 


y32 A^ 


(«) 

2/33 


/ 




yii A3 


yi2 A3 


,,('^) 

yi3 A4 


\ 


nid = 


V2\ A3 


y22 A^ 


y23 A^ 






\ ^31 A3 


y32 A2 


,,('^) 've 
y33 A 


/ 



(41) 



Vd 



(42) 



where and y^^^ (z,j = 1,2,3) are the sum of all the different terms appearing in the 
superpotential, all of them are expect to be of order one. We note that the contribution 
of d^Qj,[^x)\i^^d vanishes with the LO vacuum alignment, accordingly the (13) element of 
the down quark mass matrix arise at order 1/A^. Diagonalizing the above quark mass 
matrices in Eq. fHT]) and Eq. fll2]) with the standard perturbation technique, we obtain the 
quark masses as follows 



rrir. 



nit — 

nid ~ 

rris ^ 
rrih — 



2/11 

2/22'^^. 

(n) 
2/33 ^« 

{d)vlv^ 

yli^vd 

7)2 

2/22 ^^<i 

(d) ve 
2/33 ^Wrf 



(43) 
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We see that the quark mass hierarchies are correctly produced if the VEVs w^, v^, vq and i^a 
are of order (!?(A^A) as well. This is consistent with our naive expectation that all the VEVs 
should be of the same order of magnitude. Note that the quark mass hierarchies are generated 
through the spontaneous breaking of the flavor symmetry instead of the FN mechanism. It 
is obvious that the mass hierarchies between top and bottom quark mainly come from the 
symmetry breaking parameter vg/A, and tan/5 = ^ should be of order one in our model. 
Comparing with the tau lepton mass rrir predicted in Eq. (fT8i) . we see that rrir and rrih are of 
the same order, this is consistent with the b—r unification predicted in many grand unification 
models. 

For the quark mixing, the CKM matrix elements are estimated as 

Vud ^ Vcs ^Vtb^l 

(u) {d) (m) (d) „,3 
T.* _ 2/22 1/31 - ^211/32 



Y* ^ T/ ^32 



'td — 



(44) 



We see that the correct orders of the CKM matrix elements are reproduced with the exception 
of Cabibbo angle. The Kts(or Vcd) is the combination of two independent contributions of 

order A^, we need an accidental enhancement of the combination I ^ j^j of order 1/Ac 

in order to obtain the correct Cabibbo angle. 



5 Phenomenological implications 

In the following we shall study the constraints on the model imposed by the observed values 
of Am^gj = vn^ — ml and Am^^^ = \m'^ — m^(m2)|. The important physical consequences 
of our model are investigated in details, and the corresponding predictions are presented. In 
this section we mainly concentrate on the neutrino sector. We assume that the right handed 
neutrino mass M is much smaller than the cutoff scale A of the theory, then the light neutrino 
masses are dominantly generated via the See-Saw mechanism. 



5.1 The neutrino mass spectrum 

According to Eq. fl28l) . the light neutrino mass spectrum is controlled by two parameters 
a and b, which are in general both complex numbers. For convenience, we define 

- = i?e^* (45) 
a 

with R = I ^ I . As we will see in the following, all the low energy observables can be expressed 
in terms of only three independent quantities: the ratio R, the relative phase $ between a 



13 




Normal hierachy 



0.15 0.20 



Figure 1: The variation of cos<f> with respect to the lightest neutrino mass mi for the normal hierarchy 
spectrum. In Fig. [TJi, cos<i> is the taken to be the expression in Eq. ((48|) . and Fig. [T}d corresponds to the value 
of cos$ in Eq.ign]). 



and b, and the lightest neutrino mass. Experimentally, only two spectrum observables Am^ 



sol 



and Am^j^ have been measured, therefore the light neutrino mass spectrum can be normal 
hierarchy (NH) or inverted hierarchy(IH). The ratio between Am^^; and Am^j^ is given by 

24i?(l + 9i?2)|cos$| ^"^^^ 
mil for both the NH and IH neutrino spectrums for 



Aml^i _ 15 - 81i?^ - 18/22 _ sqr"^ cos^ $ + 12i?(l + 9/?^) cos $ 



where we have taken Am^^^ 



convenience. Moreover, we have the following relationships for the neutrino masses 
16mf 



16m| 

2 

rrio 



1 + Sli?^ + 18/2^ + 36i?2 cos^ $ - 12i?(l + 9/?^) cos $ 
1 + 81i?^ + 18i?2 + 36i?2 cos^ $ + 12i?(l + 9R^) cos $ 



(47) 



Then the parameters R and cos $ can be expressed in terms of light neutrino mass as follows 

R 



1 /2(m3+mi) 



m2 



- 1 



COS $ 



(4J 



7712 


/2(m3+mi) 




V "2 



or 



cos $ 



1 / 2|m.3-r?ii| 



"12 



- 1 



(49) 



m2|m3-mi| / 2|»Ti^-mjj _^ 



These results hold for both the normal hierarchy and inverted hierarchy spectrum. In the 
case of normal hierarchy, m2 and m^, can be expressed as functions of the lightest neutrino 
mass: m2 = y^m^~+~Am^ and = ^/mf + Am^^^. For the inverted hierarchy, is 
the lightest neutrino mass, the remaining two masses are mi = \/ m| + Am^^^ and m2 = 
\/m1 + Arn^^i + Am^j^. As a result, taking into account the experimental information on 
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Inverted hierachy 



0.10 
nil (eV) 
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Figure 2: cos$ as a function of the lightest neutrino mass for the inverted hierarchy spectrum. Fig. [5^ 
and Fig. [^h arc for the cos 4> values in Eq. PB)) and Eq. P5|) respectively. 



Am^gj and Am^^^, there is only one real parameter undetermined, and it is chose to be 
the lightest neutrino mass m/(mi or ma) in the present work. Hence our model is quite 
predictive. We display cos $ as a function of the lightest neutrino mass in Fig{T] and Figl2] 
for the normal hierarchy and inverted hierarchy respectively, where the best fit values of 
Am^g; = 7.65 X 10~^ eV^ and Am^j.^ = 2.40 x 10~^ eV^ have been used. For the solution 
of R and cos$ shown in Eq. fj49|) . we can clearly see that the corresponding value of | cos$| 
would be larger than 1 in the case of both normal hierarchy and inverted hierarchy spectrum. 
Furthermore, we have verified that | cos $| is always larger than 1 for the 3a range of Am^^j and 
Am^j^, therefore the solution in Eq.( l49i) can be disregarded thereafter. From the condition 
I cos$| < 1, we obtain the following constraints on the lightest neutrino mass, 

nil > 0.011 eV, Normal hierarchy 

ms > 0.0 eV, Inverted hierarchy (50) 

For the NH spectrum, we have a lower bound on mi, which is satisfied for $ = 0. The 
corresponding values of cos*!* are positive, then $ is in the range of ~ 7r/2 or 37r/2 ~ 27r. 
In the case of IH, cos$ is negative so that $ varies between tt /2 to 37r/2. From Fig. [2] we 
can see that cos $ is very close to -1 for tending to zero, the lightest neutrino mass ma is 
less constrained. 



5.2 Neutrinoless double beta decay 

Neutrinoless double beta decay (0i^2/?) is a sensitive probe to the scale of the neutrino 
masses, it is a very slow lepton-number-violating nuclear transition that occurs if neutrinos 
have mass and are Majorana particles. The rate of Oz/2/5 decay is determined by the nu- 
clear matrix elements and the effective Oj/2/3-decay mass |m^/}|, which is defined as |m/3/3| = 
I J2ki^PMNs)lk''^k\- In the present model it is given by 



K^| = |2a2 + 662-4a6|^ 



m2 



l + 9R^ + AR^ + QR^ cos(2$) - 12R^ cos $ - 4i? cos $ 



1/2 



(51) 
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By using Eq.( l48l) . we can express \Tnpp\ in terms of the lightest neutrino mass m/, the corre- 
sponding results are shown in Fig|31 The vertical line represents the future sensitivity of the 
KATRIN experiment [41], the horizontal ones denote the present bound from the Heidelberg- 
Moscow experiment [42] and the future sensitivity of some Qv2(3 decay experiments, which are 
15 meV, 20 meV and 90 meV, respectively of CUORE [43], Majorana [44] /GERDA III [45] 
and GERDA II experiments. From Fig. [3] we conclude that for the allowed values of m^, the 
predictions for mpp approach the future experimental sensitivity. For the NH spectrum, the 
effective mass can reach a very low value about 7.8 meV. Whereas the lower bound for 
m^^ is approximately 44.3 meV in the case of IH. A combined measurement of the effective 
mass mpp and the lightest neutrino mass can determine whether the neutrino spectrum is NH 
or IH in our model. 



1 




m(eV) 



Figure 3: as a function of the lightest neutrino mass mj, the sohd and dashed hues represent the NH 

and IH cases respectively. 



5.3 Beta decay 

One can directly search for the kinetic effect of nonzero neutrino masses in beta decay 
by modification of the Kurie plot. This search is sensitive to neutrino masses regardless of 
whether the neutrinos are Dirac or Majorana particles. For small neutrino masses, this effect 
will occur near to the end point of the electron energy spectrum and will be sensitive to the 
quantity = [^^ \ {UpMNs)ek\^'^'k^^'^ ■ For the present model, we have 

mp = ^{2ml + mlY/^ (52) 

This result holds for both the NH and IH spectrums. In FigJUwe plot mp versus the lightest 
neutrino mass m;, the horizontal line represents the future sensitivity of 0.2 eV from the 
KATRIN experiment. 
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Figure 4: Variation of mp with respect to the hghtest neutrino mass nii, the soHd and dashed hnes represent 
the NH and IH spectrum respectively. 

5.4 Sum of the neutrino masses 

The sum of the neutrino masses Ylk "^fc constrained by the cosmological observation. In 
Fig. [5] we display the sum of the neutrino masses as a function of the hghtest neutrino mass 
nil. The vertical line denotes the future sensitivity of KATRIN experiment, and the horizontal 
lines are the cosmological bounds [46]. There are typically five representative combinations 
of the cosmological data, which lead to increasingly stronger upper bounds on the sum of 
the neutrino masses. We show the two strongest ones in Fig. [51 The first one at 0.60 eV 
corresponds to the combination of the Cosmic Microwave Background (CMB) anisotropy 
data (from WMAP 5y [47], Arcminute Cosmology Bolometer Array Receiver (ACBAR) [48], 
Very Small Array (VSA) [49], Cosmic Background Imager (CBI) [50] and BOOMERANG 
[51] experiments) plus the large-scale structure (LSS) information on galaxy clustering (from 
the Luminous Red Galaxies Sloan Digital Sky Survey (SDSS) [52]) plus the Hubble Space 
Telescope (HST) plus the luminosity distance SN-Ia data of [53] and finally plus the BAO 
data from [54]. The second one at 0.19 eV corresponds to all the previous data combined to 
the small scale primordial spectrum from Lyman-alpha (Lya) forest clouds [55] . We see that 
the current cosmological information on the sum of the neutrino masses can hardly distinguish 
the NH spectrum from the IH spectrum. 

5.5 The Majorana CP violating phases 

In the standard parametrization [5], the lepton PMNS mixing matrix is defined by 

(C12C13 S12C13 5136"*'' \ 

-S12C23 - ci2S23Si3e*'' C12C23 - si2S23Si3e*'' S23C13 diag(l, 6*°''/^ 6*"''/^) 
iS iS I 

■S12S23 ~ Ci2C23Si3e — C12S23 — Si2C23Si3e C23C13 / 

(53) 

where Cy = cos 9ij, Sij = sin^ with 9ij G [0, 7r/2], 5 is the Dirac CP violating phase, 021 and 
031 are the two Majorana CP violating phases, all the three CP violating phases 6, 021 and 
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Figure 5: The sum of the neutrino masses mk versus the hghtest neutrino mass m;, the sohd and dashed 
hues represent the NH and IH spectrum respectively. 



asi are allowed to vary in the range of ~ 27r. Recalling that the leptonic mixing matrix is 
given by Eq. fE^ at LO, in the standard parametrization it is, 

UPMNS = e-*"^/' diag(l, 1, -1)Utb diag(l, e*("i-"2)/2^ ei(-i-"3)/2) ^54^ 

where the overall phase e"*"^/^ can be absorbed into the charged lepton fields. Comparing 
Eq. flS^ with Eq. fl53l) . we can identify the two CP violating phases as 

"21 = «! — 02, asi = «! — as (55) 

Similar to other low energy observables, 0121 and a^i can be written as functions of the 
parameters R and $, 



(56) 



(57) 



_ 9-R2 sin(2't>)-6/?siii<t> 
Sma2l — i+9R2_6i{cos* ' 

_ l+9fi2 cos(2't')-6i^ cos $ 
COS 021 - l+gija-eRcos* 

_ 12-R(l-9fl2)sin<t> 
Smasi — i+81/?4+18R2-36fl2 cos2 * 

_ l+81ii!"^-36fi2+i8i;2 cos(2't') 
COSaai — i+81R4+i8i?2_36ij2 cos2 $ 

Note that the relations between R, $ and the light neutrino masses are displayed in Eq. lHHl) . 
In contrast to other low energy observables such as the light neutrino masses, m^^ and 
etc., the Majorana phases 021 and 031 depend on both cos $ and sin $, and not only on cos $. 
In Figj6]we show the behavior of the Majorana phases 021 and 031 with respect to the lightest 
neutrino mass mi, where we choose sin$ > for illustration. 

It is well-known that the See-Saw mechanism provides an elegant explanation for the 
smallness of the neutrino mass, meanwhile the baryon asymmetry may be produced through 
the out of equilibrium CP violating decays of the right handed neutrino u'^. As is shown 
in Eq. (l23|) . at LO the heavy neutrino masses are exactly degenerate in our model, then 
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leptogenesis can be naturally implemented via the so-called resonant leptogenesis mechanism 
[56]. Since more other subtle issues are involved in the resonant leptogenesis, the analysis of 
whether the observed baryon asymmetry can be naturally generated in our model is beyond 
the range of the present paper, which will be discussed in future work [57]. 




Figure 6: The dependence of the Majorana CP violatmg phases a2i and a^i on the hghtest neutrmo mass 
m;. Sohd and dashed hnes refer to 0121 and a^i respectively. Fig[5^ corresponds to the NH mass spectrum, 
and Figl6)D is for the IH case, where sin $ is taken to be positive. 



The phenomenological consequences for the LO order predictions of our model have been 
analyzed . As we shall show in the next section that our model gets corrections when higher 
dimensional operators are included in the Lagrangian. These corrections modify the leading 
order predictions by terms of relative order A^, hence the results presented so far are still 
correct approximately. However, we would like to note that due to the next to the leading order 
contributions, a cancallation could be present for the effective Oi/2/3-decay mass, consequently 
mpp could reach zero in the NH case. 

6 Next to the leading order corrections 

The results of the previous section hold to first approximation. At the next to leading 
order (NLO), the superpotentials we, w^,, Wu and Wd are corrected by higher dimensional 
operators compatible with the symmetry of the model, whose contributions are suppressed by 
at least one additional power of A. The residual Klein four and Z3 symmetry in the neutrino 
and the charged lepton sectors at LO would be broken completely by the NLO contributions. 
The NLO terms in the driving superpotential leads to small deviation from the LO vacuum 
alignment. The masses and mixing matrices are corrected by both the shift of the vacuum 
configuration and the NLO operators in the Yukawa superpotentials W£, w^, and Wd- In 
the following, the NLO corrections to the vacuum alignment and the mass matrices will be 
discussed one by one, and the resulting physical effects are studied. 

6.1 Corrections to the vacuum alignment 

The NLO operators of the driving superpotential Wy and the corresponding corrections to 
the LO vacuum alignment in Eq.Q and Eq. (fT2!) are discussed in the Appendix B in details. 
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The inclusion of the higher dimensional operators results in the shift of the VEVs of the fiavon 
fields, the vacuum configuration is modified into 

(0) =ve + 5ve (58) 

where and are still undetermined, and the VEV f a is not corrected by the NLO terms. 
All the corrections are suppressed by 1/A, and the shift of (0) turns out to be proportional 
to its LO VEV. Since all the VEVs are required to be of order O(A^A), we expect these 
corrections would modify the LO VEVs by terms of relative order A^. 



6.2 Corrections to the mass matrices 

The corrections to the fermion mass matrices originate from two sources: the first is the 
higher dimensional operators in the Yukawa superpotentials we, Wu, Wu and Wd, and the second 
is the deviation from the LO vacuum alignment, which is induced by the NLO terms in the 
driving potential. As a result, at NLO the mass matrices are the sum of the contributions of 
higher dimension operators evaluated with the insertion of the the LO VEV, and those from 
the LO superpotentials evaluated with the NLO VEVs. 

For the charged leptons, the superpotential we is corrected by the following sisteen NLO 
operators 

e^(V)i.^/^d, eV^\)i,ehd, e\i^x\,Ohd. e\lx\,Ohd, 
e'{e^r]^)i,Ahd, e'{ixV^)i,Ahd, e''{i^7]^)i,Ahd, e'{ixV^)i,Ahd, 
e^(£<^</.2)i,A/irf, e%ex(p')i,Ahd, /i^(V)u^/id, f^Vx')uOhd, 
^i\i^xKOhd. ^lV<P%Ahd. fi''{iv<l>)uAhd, r'^iixh^Ohd (59) 

Taking into account the contributions of the modified vacuum alignment at NLO, each di- 
agonal entry of the charged lepton mass matrix receives a small correction factor, while the 
off-diagonal entries become non-zero and of the order of the diagonal term in each row mul- 
tiplied by e, which parameterizes the ratio VEV/ A with order 0{Xl). Then we have 

m^^e^ m22£ ""^23^^ (60) 
mii£ ml^e m^^ / 

where the coefficients mlJi,j =1,2,3) are order one unspecified constants. The hermitian 
matrix m^rrii is diagonalized by the unitary matrix Ui, which exactly corresponds to the 
transformation of the charged leptons used to diagonalize m^, 

UjmlmcUi ~ diag(|m^^e^p, |m22£:^P, \'ml^e\'^)vl (61) 

The charged lepton masses are modified by terms of relative order e with respect to LO results, 
consequently the NLO corrections don't spoil the charged lepton mass hierarchies predicted 
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at LO. The unitary matrix Ui is approximately given by 



( 



\ 



'^22 

'^31 c 
i ^ 

■'-as 



'"■22 

1 

"^33 



'^32, 



(62) 



Then we turn to the neutrino sector. The NLO correction to the Majorana masses of the 
right handed neutrino arises at order 1/A, the corresponding higher dimensional operator is 
{i''^u'^)i-^6'^, whose contribution can be completely absorbed into the redefinition of the mass 
parameter M. The NLO corrections to the neutrino Dirac couplings are 



A 



yu2 
A 



uOhu 



(63) 



where Sr] and 50 represent the shifted VEVs of the fiavons rj and respectively. Through 



redefining the LO parameters a 



a — x^i^^jT- and b 



are 



5m 



D 




5i 



6H 

S2 



yu2- 



A ' 



^2 





the NLO corrections to 



(64) 



where 5i = Uui—^ + '^Xi^i^j^ and 62 = Xi^2-j^- We notice that both 61 and 62 are of order 
ea (or eb). Therefore, the NLO corrections to the light neutrino mass matrix are given by 



K 
M 



-2(a - b)6i 
-(2a + b)5i - 66^2 
-6(5i - 662) 



N 

(2a + b)6i - 6M2 
2b{6i + 3^2) 
-(2a + b)6i 



-b{6i - 6S2) 
-(2a + b)Si 
-2(a - b)5i - Qb62 



(65) 



Diagonalizing the modified light neutrino mass matrix, we obtain the neutrino masses to LO 
in (5i.2 as follows. 



nil 



m2 



a - 3bf + (a - 3b)5i 
,,2 



\M\ 



a6i 



[a + 3bf 



\M\ 



3b)Si 



\M\ 



(66) 



The PMNS matrix becomes Upmns = UjUu, where Ui associated with the diagonalization of 
the charged lepton mass matrix is given by Eq. fl62l) . and the unitary matrix f/j, diagonalizes 
the neutrino mass matrix m,^ + 6mi, including the NLO contributions. The parameters of the 
lepton mixing matrix are modified as 



Ue3\ = 



y/2 6(|a|2 + 9|6|2)(a6* + a*6) 



[(a + 3bYia*5l + 6b*5; 



3b*y{a6i + 6M2)] 



+ 



m 



21 



m 



22 



31 . 



m' 



33 



21 



sin^^ 



12 



sin^ e~ 



23 
1 

+ 2 



1 

3 L 

1 

2 



m 



21 



31 



e — 



m 



22 



33 



m 



21 



£ - 



m 



31 



22 ' ^ ""^33 

ab{a*6l + 66*5;) + a*b*{a6i + 6M2) 



2(|a|2 + 9|6|2)(a6* + a*6) 



m 



32 



33 



32 , 



33 



(67) 



We see that the neutrino masses and mixing angles receive corrections of order with 
respect to LO results. The value of sin^6'i2 is still within the 3a range of global data fit, 
and the corrections to both 6*23 and ^13 are within the current data uncertainties as well. In 
particular, a non- vanishing 613 of order is close to the reach of the next generation neutrino 
oscillation experiments and will provide a valuable test of model. 

The NLO corrections to the mass matrices in quark sector have been analyzed following 
the same method as that for the lepton sector. Since every entry of the mass matrices 
rriu and 171^ in Eg. (141 p and Eg. (1421) is nonvanishing, the NLO contributions leads to small 
corrections of relative order A^ in each entry. Conseguently the guark masses and mixing 
angles are corrected by terms of relative order A^ with respect to LO results, the successful 
LO predictions are not spoiled. 



7 Conclusion 

We have constructed a SUSY model for fermion masses and flavor mixings based on the 
flavor symmetry 6*4 x Z3 x Z4, the neutrino masses are assumed to be generated through the 
See-Saw mechanism. At LO the 5*4 symmetry is broken down to Klein four and Z3 symmetry 
in the neutrino and charged lepton sector respectively, this breaking chain exactly leads to 
the TB mixing. It is remarkable that the mass hierarchies among the charged leptons are 
controlled by the spontaneous breaking of the flavor symmetry. We further extend the flavor 
symmetry to the guark sector, where the 5*4 symmetry is completely broken. The correct 
orders of guark masses and CKM matrix elements are generated with the exception of the 
mixing angle between the first two generations, which reguires a samll accidental enhancement. 

We have carefully analyzed the NLO contributions due to higher dimensional operators 
which modify both the Yukawa couplings and the LO vacuum alignment, and we have verified 
that all the fermino masses and mixing angles are corrected by terms of relative order A^ 
with respect to the LO results. As a result, the successful LO predictions are not spoiled. 
Particularly we expect the mixing angle ^13 would be of order A^, it is within the sensitivity 
of the experiments which are now in preparation and will take data in the near future [58,59]. 
Precise measurement of ^13 is an important test to our model. 

The phenomenological conseguences of our model are analyzed in detail. The low energy 
observables including the neutrino mass sguared difference, neutrinoless double decay, beta 
decay, the sum of the neutrino masses and the Majorana CP violating phases are considered. 
All the low energy observables can be expressed in terms of three independent parameters: 
the ratio R = \b/a\, the relative phase $ between a and b and the lightest neutrino mass 
mi. Once the parameters are fixed to match Am^^j and Am^^^, there is only one parameter 
left, which is chose to be mi in the present work. Both the normal and inverted hierarchy 
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neutrino spectrum are allowed in our model. For normal hierarchy there is a low bound on 
nil of approximately 0.011 eV. In the case of inverted hierarchy, m^, is less constrained and 
we only obtain the trivial constraint that should be positive. The lower bounds of the 
effective mass are approximately 7.8 meV and 44.3 meV for the NH and IH spectrum 
respectively. A combined measurement of m^^ and the lightest neutrino mass can distinguish 
the NH from the IH spectrum. The Majorana CP violating phases depend both on cos $ and 
sin$, whereas only cos$ is involved in other low energy observables. It is remarkable that 
the right handed neutrino masses are exactly degenerate at LO, the baryon asymmetry may 
be generated via the rosonant leptogenesis. 
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Appendix A: Representation matrices of the S4 group 



In this appendix, we explicitly show the representation matrices of the S4 group for the 
five irreducible representations. The matrices for the generators S and T depend on the 
representations as follows 

li, 
I2, 



3i, S 



S 




T 




where cu = e^'^^l'^ = (—1 + v^)/2. In the identity representation li, all the elements are 
mapped onto the number 1. In the antisymmetric representation I2, the group elements 
correspond to 1 or -1 respectively for even permutation and odd permutation. For the 2 
representation, the representation matrices are as follows 



1 
1 



C2 : STS' = T^S=[ ° V TSTS^ = TST=(^^ M , ST' = S'TS =(^20 



1 
1 



C: T = S'T = TS' = S'TS' =i M = S'T' = STS = T'S' ^ ° 



4.1-01-10 -010 - \ 2/'- -- - - in 

C5: S = S'=(^^ l), T'ST = TS=( ° V ST = TST' ^ ° 



10/' V cu^ / ' V 



For the 3i representation, the representation matrices are 
V 1 / 

2 

S'TS = I cuM , TST = ^ I 2 2 -1 I , = ^ I 2cu" 2 -co 

2u -co' 2 
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/ 


T^S^ = 


ll 

"A 
\ 


' -1 


2uj 


2u^ 


2uj 




2 




2 


— UJ 


(;\ 


2 


2 


2u^ 






V 2uj 


2u 


— UJ 




2 


2 


2uj^ 


2u^ 




^ 2uj 


—UJ 


2uj 




2uj^ 


2uj 




2uj^ 


— UJ 






2uj 




2uj 
2 

-UJ^ 



,ST 




Since the signs of the generator S are opposite in 3i and 82 representations, the represention 
matrices for the 82 representation can be found from those of the 81 representation: the 
matrices are exactly the same for Ci, C3 and C4 classes, whereas they are the opposite for C2 
and C5. 



Appendix B: NLO corrections to the vacuum alignment 

In this appendix, we will analyze the NLO corrections to the vacuum alignment induced 
by the higher dimensional operators. At NLO the driving superpotential dependent on the 
driving fields is modified into 

+ Aw^ (68) 

where is the LO contributions shown in Eq.([7]), which are of dimension three. Aw^ is most 
general set of terms suppressed by one power of the cutoff, which are linear in the driving 
fields and are invariant under the symmetry of the model. Concretely Aw^ is given by 

A^. = ^ [ p^of + ^2 + E + E + E '^^t] (69) 

i=l i=l 1=1 1=1 i=l 

where pi, Xj, Cj and Si are order one coefficients, [O^Of .OlOl.Of} are the complete 
set of invariant operators of dimension four, 

Ot = (^°(^x)3.)i. 01 = (^°(#)3ji, A (70) 
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0[ = d\r,\,, 0', = 9\4>'),„ Ol = d\^<p\, (72) 
O'r = mv')2)iA Ol = {v'{<P%)i.e (73) 

Ot = (0°(0^)3ji.^, Ot = i<P%V<PK)uO (74) 

The NLO superpotential Awy results in shift of the LO VEVs, then the vacuum configuration 
is modified into 

{r]) = {v.^, Vr, + SVri^), {(p) = {v^ + Sv^^,v^ + Sv^^,v^ + 5v^^), 

(9) =ve + Svg (75) 

Note that the VEV doesn't receive correction at NLO. Similar to section [3l the new 
vacuum configuration is obtained by imposing the vanish of the first derivative of Wy + Awy 
with respect to the driving fields ip^, C,'^, 6^, rf and (fP. Only terms linear in the shift 5v are 
kept, and terms of order 6v/K are neglected, then the minimization equations become, 

i2giv^ + gsVx)6v^.^ + {2g2V^ - gsv^)6v^^ = 

2g,6v^,+p,^ = 

{2giv^ - g2.v^)5v^^ + {2g2V^ + gzv^)5v^^ = 

gA{v^6v^^ + v^6v^.^) + —{2x3Vav'^ + Sx^vav^ + X5vD = 

nve6ve + + ^hvy^) = (76) 

Solving the above linear equations, we obtain 

93Vx - 2giVp 



Pi Vx^e 
'2gi A 

2g2Vx - 93Vp 



= 2,4A(2,,.^ + ,3.,.,) ['"^"^"-' + '"^"^"^ + 

= -2^;JSrS;;^^'"^^"^' + '"^^^^^ + "^^^^ 

^ve = -r^ -r-^ (77) 

K Ave K Af e 

where Sv^^.s^ ^"^xs ^'^s are of order 1/A, 6vxi is undetermined, and it is expected to be 
suppressed by 1/A as well. The minimization conditions for the shift 6v^^ and 6v^-^^ ^ 3 are 

2f2V4>{5v^^ + 5^02 + 5v^^) + ei^^ + 3e2^^ = 

9 2 

2fiVr,5vr,2 + 2f2V^{6v^, + 6v^^ + 6v^^) - ei— 3e2 



A ^ A 

fsiv^Sv^^ - v^5v^^ - v^6vn^) + 2^2 ^''^^^ = 

hiVr^dv^^ - Vr,6v^2 - V4>5Vr^^) + 2^2 ^''^^^ = 

h{vr^5v^^ - Vr,6v^^ - v<p6v^^) + 2^2 ^''^^^ = (7J 
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The solutions to the above equations are 



2S2 VnVg 



h A 



X - X K A ^1 ^^^^ ^2 v^ve , . 

Sv, = Sv,, = = ^3 = -g^^ - (79) 

We notice that (0) acquires 0{1/A) corrections in the same directions, and the shift Svri2 is 
in general non-zero. Since all the VEVs approximately are of the same order O(A^A) at LO, 
we expect ^y^y ~ ^^^^ ~ A^. Therefore the LO vacuum alignment in Eq.Q and Eq. f[T21) is 
stable under the NLO corrections. 



27 



References 



[1] T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 10, 113011 (2008) 
|arXiv:0808.2"0T6l [hep-ph]]; M. Maltoni and T. Schwetz, larXiv:0812.3161l [hep-ph]. 



[2] M. C. Gonzalez-Garcia and M. Maltoni, Phys. Kept. 460, 1 (2008) [arXiv: 0704 .18001 
[hep-ph]]. 

[3] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002), 



hep-ph/0202074- P. F. Harrison and W. G. Scott, Phys. Lett. B 535, 163 (2002), 



hep-ph/0203209: Z. Z. Xing, Phys. Lett. B 533, 85 (2002), hep-ph/0204049, X. G. He 



and A. Zee, Phys. Lett. B 560, 87 (2003), hep-ph/0301092, 



[4] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. 
Theor. Phys. 49, 652 (1973). 

[5] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 

[6] S. F. King, JHEP 0508 (2005) 105 |arXiv:hep-ph/0506297j; I. de Medeiros Varzielas and 



G. G. Ross, arXiv:hep-ph/0507176; I. de Medeiros Varzielas, S. F. King and G. G. Ross, 
Phys. Lett. B 644 (2007) 153 |arXiv:hep-ph/0512313] ; I. de Medeiros Varzielas, S. F. King 
and G. G. Ross, Phys. Lett. B 648 (2007) 201 |arXiv:hep-ph/0607045l ; S. F. King and 
M. Malinsky, JHEP 0611 (2006) 071 |arXiv:hep-ph/0608021| ; F. Feruglio and Y. Lin, 
Nucl. Phys. B 800, 77 (2008) iarXiv:0712. 15281 [hep-ph]]. . 



[7] M. Leurer, Y. Nir and N. Seiberg, Nucl. Phys. B 420, 468 (1994) |arXiv:hep-ph/9310320 



M. Leurer, Y. Nir and N. Seiberg, Nucl. Phys. B 398, 319 (1993) |arXiv:hep-ph/9212278| ; 



L. E. Ibanez and G. G. Ross, Phys. Lett. B 332, 100 (1994), arXiv:hep-ph/9403338 P. Bi- 



netruy and P. Ramond, Phys. Lett. B 350, 49 (1995), arXiv:hep-ph/9412385 : E. Dudas, 



S. Pokorski and C. A. Savoy, Phys. Lett. B 356, 45 (1995), |arXiv:hep-ph/9 504292| 

G. G. Ross, "Models of fermions masses". Prepared for Theoretical Advanced Study In- 
stitute in Elementary Particle Physics (TASI 2000): Flavor Physics for the Millennium, 
Boulder, Colorado, 4-30 Jun 2000; G. AhareUi. [arXiv:0711.016T1 [hep-ph]. 



[9] E. Ma and G. Rajasekaran, Phys. Rev. D 64 (2001) 113012, 'arXiv:hep-ph/0106291 



E. Ma, Mod. Phys. Lett. A 17 (2002) 627 [arXiv:hep-ph/0203238] ; K. S. Babu, E. Ma and 
J. W. F. Valle, Phys. Lett. B 552 (2003) 207, ]arXiv:hep-ph/0206292l M. Hirsch, J. C. Ro- 



mao, S. Skadhauge, J. W. F. Valle and A. Villanova del Moral, arXiv: hep-ph/ 03 12 244 



Phys. Rev. D 69 (2004) 093006 |arXiv:hep-ph/0312265| ; E. Ma, Phys. Rev. D 70 
(2004) 031901; Phys. Rev. D 70 (2004) 031901 |arXiv:hep-ph/0404199| ; New J. Phys. 
6 (2004) 104 |arXiv:hep-ph/0 405152]; 'arXi v:hep-ph/0409075t S. L. Chen, M. Frige- 
rio and E. Ma, Nucl. Phys. B 724 (2005) 423 |arXiv:hep-ph/ 0504181|; E. Ma, Phys. 
Rev. D 72 (2005) 037301 |arXiv:hep-ph /0505209l ; M. Hirsch, A. Villanova del Moral, 
J. W. F. Valle and E. Ma, Phys. Rev. D 72 (2005) 091301 [Erratum-ibid. D 72 (2005) 
119904] |arXiv:hep-ph/0507148| . K. S. Babu and X. G. He, |arXiv:hep-ph/0507217; 
E. Ma, Mod. Phys. Lett. A 20 (2005) 2601 [arXiv:hep-ph/0508099|; A. Zee, Phys. 



28 



Lett. B 630 (2005) 58 |arXiv:hep-ph/0508278l ; E. Ma, Phys. Rev. D 73 (2006) 057304 
|arXiv:hep-ph/0511133] ; X. G. He, Y. Y. Keum and R. R. Volkas, JHEP 0604 (2006) 



039 |arXiv:hep-ph/0601001| ; B. Adhikary, B. Brahmachari, A. Ghosal, E. Ma and 



M. K. Parida, Phys. Lett. B 638 (2006) 345 |arXiv:hep-ph/060305 9|; E. Ma, Mod. 
Phys. Lett. A 21 (2006) 2931 |arXiv:hep-ph/0607190| ; Mod. Phys. Lett. A 22 (2007) 



101 [arXiv:hep-ph/0610342| ; L. Lavoura and H. Kuhbock, Mod. Phys. Lett. A 22 (2007) 



181 [arXiv:hep-ph/06 10050]; E. Ma, H. Sawanaka and M. Tanimoto, Phys. Lett. B 641, 



301 (2006), hep-ph/0606103i S. F. King and M. Mahnsky, Phys. Lett. B 645 (2007) 351 



|arXiv:hep-ph/06 1 0250] ; S. Morisi, M. Picariello and E. Torrente-Lujan, Phys. Rev. D 
75 (2007) 075015 |arXiv:hep-ph/0702034 |; M. Hirsch, A. S. Joshipura, S. Kaneko and 
J. W. F. Valle, Phys. Rev. Lett. 99, 151802 (2007) [arXiv:hep-ph/0703046 |. F. Yin, 
Phys. Rev. D 75 (2007) 073010 [arXiv: 0704. 38271 [hep-ph]]: F. Bazzocchi, S. Kaneko and 
S. Morisi, JHEP 0803 (2008) 063 |arXiv:0707.3032l [hep-ph]]. F. Bazzocchi, S. Morisi and 
M. Picariello, Phys. Lett. B 659 (2008) 628 |arXiv:0710.2928l [hep-ph]]; M. Honda and 
M. Tanimoto, Prog. Theor. Phys. 119 (2008) 583 [arXiv:0801.0181 [hep-ph]]; B. Brah- 
machari, S. Choubey and M. Mitra, Phys. Rev. D 77 (2008) 073008 [Erratum-ibid. D 
77 (2008) 119901] [arXiv:0801 .3554' [hep-ph]]: F. Bazzocchi, S. Morisi, M. Picariello and 
E. Torrente-Lujan, J. Phys. G 36 (2009) 015002 [arXiv:0802. 1693 1 [hep-ph]] : B. Adhikary 
and A. Ghosal, Phys. Rev. D 78 (2008) 073007 jarXiv:08 03.3582 [hep-ph]]; M. Hirsch, 
S. Morisi and J. W. F. Valle, Phys. Rev. D 78 (2008) 093007 ja?Xrv:0804. 15211 [hep-ph]]: 
Y. Lin, Nucl. Phys. B 813, 91 (2009) [arXiv: 0804. 28671 [hep-ph]]; P. H. Frampton and 
S. Matsuzaki. [arXiv:0806.4592l [hep-ph]: F. Feruglio, C. Hagedorn, Y. Lin and L. Merlo, 
Nucl. Phys. B 809, 218 (2009) |arXiv:0807.3T60l [hep-ph]]; F. Bazzocchi, M. Frigerio 
and S. Morisi, Phys. Rev. D 78, 116018 (2008) [arXiv:0809.3573' [hep-ph]]; S. Morisi, 
Phys. Rev. D 79, 033008 (2009) [arXiv:0901.1080 [hep-ph]]; P. Ciafaloni, M. Picariello, 
E. Torrente-Lujan and A. Urbano, Phys. Rev. D 79, 116010 (2009) [arXiv:0901.2 236 [hep- 
ph]]; M. C. Chen and S. F. King, JHEP 09 06, 072 (2009) |arXiv:0903.0T25l [hep-ph]]; 
T. J. Burrows and S. F. King. larXiv:0909. 14331 [hep-ph]. . 



[10] G. AltareUi and F. Feruglio, Nucl. Phys. B 720, 64 (2005), >hep-ph/0504165 



[11] G. Altarelh and F. Feruglio, Nucl. Phys. B 741, 215 (2006), hep-ph/0512103 



[12] G. Aharelh, F. Feruglio and C. Hagedorn, JHEP 0803, 052 (2008) |arXiv:0802.0"090l 
[hep-ph]]. 

[13] G. C. Branco, R. Gonzalez Felipe, M. N. Rebelo and H. Serodio, Phys. Rev. D 79, 093008 
(2009) |arXiv:0904.3076l [hep-ph]]. 

[14] G. Altarelh and D. Meloni, J. Phys. G 36, 085005 (2009) |arXiv:0905.0"620l [hep-ph]]. 

[15] E. Bertuzzo, P. Di Bari, F. Feruglio and E. Nardi. [arXiv:0908.0161l [hep-ph]. 

[16] C. Hagedorn, E. Molinaro and S. T. Petcov, JHEP 0909, 115 (2009) [arXiv:0908.0"240l 
[hep-ph]]. 



29 



[17] D. Aristizabal Sierra, F. Bazzocchi, I. de Medeiros Varzielas, L. Merlo and S. Morisi, 



[18 
[19 
[20 
[21 

[22 

[23 
[24 
[25 

[26 

[27 
[28 
[29; 
[30 
[31 

[32 
[33 



,arXiv:0908. 0907 [hep-ph]. 

R. G. Felipe and H. Serodio. [a?Xiv:0908.2947l [hep-ph]. 

G. J. Ding, Phys. Rev. D 78, 036011 (2008) |arXiv:0803.2278l [hep-ph]]. 



P. D. Carr and P. H. Frampton, arXiv:hep-ph/0701034 



F. Ferugho, C. Hagedorn, Y. Lin and L. Merlo, Nucl. Phys. B 775, 120 (2007), 



hep-ph/0702194 



M. G. Ghen and K. T. Mahanthappa, Phys. Lett. B 652, 34 (2007), larXiv:0705.0"7T4l 
[hep-ph] . 

P. H. Frampton and T. W. Kephart, JHEP 0709, 110 f2007). [ailGv:0706. 11861 [hep-ph]. 

A. Aranda, Phys. Rev. D 76, 111301 r2007). Ia^v:0707.3661l [hep-ph]. 

P. H. Frampton, T. W. Kephart and S. Matsuzaki, Phys. Rev. D 78, 073004 (2008) 
|arXiv:0807.47T3] [hep-ph]]. 

D. A. Eby, P. H. Frampton and S. Matsuzaki, Phys. Lett. B 671, 386 (2009) 
|arXiv:0810.4899l [hep-ph]]. 

E. Ma, Phys. Lett. B 632, 352 (2006) [arXiv:hep-ph/0508231| . 

F. Bazzocchi and S. Morisi. larXiv:0811. 03451 [hep-ph]. 
F. Bazzocchi, L. Merlo and S. Morisi. [arXiv:0901.2086l [hep-ph]. 

F. Bazzocchi, L. Merlo and S. Morisi. [arXiv:0902.2849 [hep-ph]. 

H. Ishimori, Y. Shimizu and M. Tanimoto, Prog. Theor. Phys. 121, 769 (2009) 
|arXiv:0812.503T] [hep-ph]]. 

G. Altarelh, F. Feruglio and L. Merlo, JHEP 0905, 020 (2009) [arXiv: 0903. 19401 [hep-ph]] . 

S. Pakvasa and H. Sugawara, Phys. Lett. B 82, 105 (1979); T. Brown, N. Deshpande, 
S. Pakvasa and H. Sugawara, Phys. Lett. B 141, 95 (1984); Y. Yamanaka, H. Sug- 
awara and S. Pakvasa, Phys. Rev. D 25, 1895 (1982) [Erratum-ibid. D 29, 2135 (1984)]; 
T. Brown, S. Pakvasa, H. Sugawara and Y. Yamanaka, Phys. Rev. D 30, 255 (1984). 

[34] D. G. Lee and R. N. Mohapatra, Phys. Lett. B 329, 463 (1994) [arXiv:hep- ph/940 3201]; 
G. Hagedorn, M. Lindner and R. N. Mohapatra, JHEP 0606, 042 (2006) 
|arXiv:hep-ph/0602244] ; Y. Gai and H. B. Yu, Phys. Rev. D 74, 115005 
(2006) |arXiv:hep-ph7060802 2|; H. Zhang, Phys. Lett. B 655, 132 (2007) 
|arXiv:hep-ph/0612214| ; Y. Koide, JHEP 0708, 086 (2007) |arXiv:0705.2275l [hep-ph]]: 
M. K. Parida, Phys. Rev. D 78, 053004 (2008) |arXiv:0804.4"57T] [hep-ph]]. 

[35] G. S. Lam, Phys. Rev. Lett. 101, 121602 (2008) [arXiv: 0804. 26221 [hep-ph]]. 



30 



[36] C. S. Lam, Phys. Rev. D 78, 073015 (2008) |arXiv:0809. 11851 [hep-ph]]. 
[37] C. S. Lam. larXiv:0907.2206l [hep-ph]. 

[38] W. Grimus, L. Lavoura and P. O. Ludl. larXiv:0906.2689l [hep-ph]. 

[39] P. Minkowski, Phys. Lett. B67 (1977) 421; T. Yanagida, in Proceedings of the Workshop 
on Unified Theory and the Baryon Number of the Universe, edited by O. Sawada and A. 
Sugamoto (KEK, Tsukuba, 1979); M. Gell-Mann, P. Ramond and R. Slansky, in Super- 
gravity., edited by P. van Nieuwenhuizen and D. Freedman (North Holland, Amsterdam, 
1979); S. L. Glashow, in Quarks and Leptons, edited by M. Levy et al. (Plenum, New 
York, 1980); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 



[40 
[41 

[42 
[43 
[44; 
[45 
[46 

[47 
[48 
[49 

[50 

[51 

[52 

[53 

[54; 



C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147, 277 (1979). 



A. Osipowicz et al. [KATRIN Collaboration], arXiv:hep-ex/0109033; see also: 



http: / /www- ik.fzk.de/ ~katrin/index.html 



L. Baudis et al, Phys. Rev. Lett. 83, 41 (1999) |arXiv:hep-ex/9902014| . 

A. Giuhani [CUORE Collaboration], J. Phys. Conf. Ser. 120 (2008) 052051. 

Majorana Collaboration, larXiv: 08 1 1 . 24461 [nucl-ex]. 

A. A. Smolnikov and f. t. G. Collaboration. larXiv:0812.4194l [nucl-ex]. 

G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, larXiv:0809.2936l [hep- 
ph]. 

WMAP Collaboration, E. Komatsu et al. larXiv:0803. 05471 [astro-ph]. 

ACBAR Collaboration, C. L. Reichardt et a/.. larXiv:0801.1491l [astro-ph]. 

VSA Collaboration, C. Dickinson et al, Mon. Not. Roy. Astron. Soc. 353, 732 (2004) 
|arXiv:astro-ph/0402498] . 

CBI Collaboration, A. C. S. Readhead et al, Astrophys. J. 609, 498 (2004) 
|arXiv:astro-ph/0402359] . 

BOOMERANG Collaboration, C. J. MacTavish et al, Astrophys. J. 647, 799 (2006) 
|arXiv:astro-ph/0507503] . 

SDSS Collaboration, M. Tegmark et al, Phys. Rev. D 74 (2006) 123507 



|arXiv:astro-ph/0608632] . 

SNLS Collaboration, P. Astier et al Astron. Astrophys. 447, 31 (2006) 
|arXiv:astro^/0510447] . 

SDSS Collaboration, D. J. Eisenstein et al, Astrophys. J. 633, 560 (2005) 
|arXiv:astro-ph/050117l] . 



31 



[55] P. McDonald et ai, Astrophys. J. Suppl. 163, 80 (2006); P. McDonald et ai, Astrophys. 
J. 635, 761 (2005). 

[56] A. Pilaftsis and T. E. J. Underwood, Phys. Rev. D 72, 113001 (2005) 
|arXiv:hep-ph/0506107] . 

[57] Gui-Jun Ding, work in progress. 



[58] F. Ardellier et al. [Double Cliooz Collaboration], arXiv:hep-ex/0606025 



Y. f. Wang, arXiv:hep-ex/0610024 



32 



